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The relativity of Global Positioning System (GPS) pseudorange measurements is explored within 
the geometrical optics approximation in the curved space-time near Earth. A space-time grid for 
navigation is created by the discontinuities introduced in the electromagnetic field amplitude by 
the P-code broadcast by the GPS satellites. We compute the world function of space-time near 
Earth, and we use it to define a scalar phase function that describes the space-time grid. We use 
this scalar phase function to define the measured pseudorange, which turns out to be a two-point 
space-time scalar under generalized coordinate transformations. Though the measured pseudorange 
is an invariant, it depends on the world lines of the receiver and satellite. While two colocated 

CO ■ receivers measure two different pseudoranges to the same satellite, they obtain correct position and 

time, independent of their velocity. We relate the measured pseudorange to the geometry of space- 
time and find corrections to the conventional model of pseudorange that are on the order of the 

s ! ■ gravitational radius of the Earth. 



PACS numbers: 95.30.Sf,91.10.-v,91.10.By,95.55.Br 



Clock synchronization algorithms play a key role in applications such as communication, message encryption, and 
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navigation. Recently, there have been a number of proposed clock synchronization algorithms based on a quantum 
information approach 0,0,0iSIEIESHl- These discussions have been based mostly on non-relativistic quantum 
mechanics. In many applications, however, clock synchronization must be performed between two nodes that are in 
relative motion, such as a jet aircraft and the ground, or between two nodes that are at different gravitational poten- 
tials, such as a satellite and the ground, or even two satellites at different altitudes. In such cases, a correct treatment 
of quantum clock synchronization must include relativistic effects from the start. The well-known tension between 

O 1 quantum mechanics and relativity theory makes this a difficult task. It is clear that the concepts of measurement, 
and transformation of measurable quantities under Lorentz transformations, play a key role in the problem of clock 

bi). synchronization for both quantum and classical schemes. However, the transformation of measurable quantities has 
not been discussed even for the case of classical clock synchronization. In this article, we analyze in some detail the 
• i-H , relativity of clock synchronization in the Global Positioning System (GPS), which is based on a classical synchroniza- 
tion scheme. In particular, we describe the transformation properties of the measured quantity called pseudorange in 
the GPS [3, [T(J, UM ■ We hope that clarifying the transformation properties of measurable quantities in classical clock 
synchronization will provide some useful insight into the problem of quantum clock synchronization. 

The Global Positioning System (GPS) is a U.S. military constellation of satellites used for time keeping, and for 
navigation of land, air and sea | 9| , llCtllllll2l |. Recently, two papers have analyzed the system of space-time coordinates 
that is used in the GPS [l3l Il4|. In this paper, we address a different but related aspect: the transformation of GPS 
pseudorange measurements. A remarkable aspect of the GPS is that a receiver need not be stationary with respect 
to the Earth's surface to obtain accurate time and position. A ship, a jet aircraft, or a low-Earth-orbit satellite can 
each compute accurate time and position, even though they have different velocities. This feature of the GPS is a 
consequence of two aspects: the signal structure of the satellite broadcasts and the special type of measurement that 
a GPS receiver makes. The GPS satellite signals set up an invariant grid of 3-dimensional space-time hypersurfaces 
(light-cones). Each hypersurface is uniquely marked by the satellite that generated the hypersurface and by the 
space-time coordinates of the event of generation of the light-cone 0, 0] . By measuring the pseudorange to four 
satellites, the GPS receiver essentially determines its position by identifying the four hypersurfaces that it intersects. 
The pseudorange measurement is independent of receiver motion, up to an additive constant. For this reason, the 
measured pseudorange may be called a Lorentz pseudo-invariant. Below, we will see that the pseudorange is actually 
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a two-point scalar under generalized coordinate transformations. The transformation properties of the pseudorange 
is a key element of the GPS, yet this subject has only briefly been mentioned in the literature |l5j . 

In order to present a coherent description of the transformation properties of the pseudorange, we must deal with 
the nature of the broadcast GPS signals, the relativistic effects that impact these signals, as well as the measurement 
process itself. Therefore, the outline of this article is as follows. In section II, we present a two-receiver thought 
experiment to clarify the concept of invariance of pseudorange measurements. Section III contains a description of 
the space-time metric in the vicinity of the Earth. Section IV describes the relativistic effects on the GPS satellite 
clocks and on the observed signals, using the metric in section III. Section V discusses the nature of the GPS broadcast 
signals, which are used to set up the space-time grid. Section VI discusses the pseudorange measurement process, 
using a mechanical analogue for a GPS receiver. The transformation properties of the pseudorange are obtained in 
section VI. Section VII contains a brief discussion of navigation and time transfer using GPS signals. A correction 
is derived to the conventional flat space-time model of pseudorange. Section IX contains a summary and comments. 



II. TWO RECEIVER EXPERIMENT 



In order to make clear the concept of transformation of pseudorange measurements we offer a simple thought 
experiment. Consider two identical GPS receivers that are in relative motion. For example, one receiver is stationary 
with respect to the Earth's surface and the other is on a jet aircraft travelling at 1000 km/hour. Each GPS receiver 
carries an identical copy of the software that is used to compute receiver position and time. Assume that each receiver 
is tracking the same four GPS satellites. See Fig. ^ Assume the two receivers' world lines cross at an event M 
in space-time. Does each receiver compute the same spatial position and time for the coincident event Ml From 
a physics stand point, the GPS satellites orbit the Earth at approximately 8.37 km/s, and there are large Doppler 
frequency shifts due to satellite and receiver motion 16]. In the GPS, the actual computation of receiver position 
and time depends on the space-time coordinates of the signal emission event (at the satellite), x\, and reception 
event (at the receiver), x % . In an Earth-centered inertial (ECI) frame, these coordinates are x\ and x % , however in the 
receiver's comoving frame the emission and reception events have different coordinates, say x'£ and x'£ . The two sets of 
coordinates, (x\, x l Q ) and (x s * , x'*), for i = 0, 1, 2, 3, are related by a Lorentz transformation, which depends on receiver 
velocity with respect to the ECI frame. The actual electromagnetic field is different in each receiver's comoving frame, 
so there is different input information into each identical measuring device (receiver computer program). How can 
the two receivers compute the same spatial position and time from different input information? The key concepts are 
the space-time grid that is created by GPS satellites and the transformation properties of the measurable quantity in 
the GPS, which is called the pseudorange. These two themes are developed in the following sections. 



III. CHOICE OF METRIC 



In order to discuss relativistic effects, a reference frame or system of coordinates must be chosen. In the weak field 
limit, the metric of space- time in the vicinity of the Earth is of the form 
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where V is the Newtonian gravitational potential and (x°, x 1 , x 2 , x 3 ), are the coordinates. The frame of reference 
in Eq. (JIJ can be thought of as an Earth-centered inertial (ECI) frame. We neglect small off-diagonal terms go a , 
a = 1,2, 3, due to the rotation of the Earth. 

In general relativity, the coordinates are mathematical entities that are never directly observed. However, it is 
useful to choose the coordinates in some physically meaningful way. The coordinates x°, x 1 , x 2 , x 3 are geocentric 
coordinates, where x 3 coincides with the Earth's axis of rotation and increa sing positive values point to North. The 
Earth is modelled as an oblate spheroid with Newtonian potential given by |l8| 



V(r,9) = - 



GM 



1- ./,(£) /Mco^li 



(2) 



where, r 2 = (x 1 ) 2 + (x 2 ) 2 + (x 3 ) 2 and 9 is the polar angle measured from the x 3 axis. In Eq. (J2J, G is Newton's 
gravitational constant, M is the mass of the earth, P%(x) — (3x 2 — l)/2 is the second Legendre polynomial, R is the 
Earth's equatorial radius, and J2 is the Earth's quadrupole moment, whose value is approximately J2 = 1.08 x 10~ 3 , 
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FIG. 1: The world lines of two receivers, 0\ and O2, are shown. The receivers are tracking the same four GPS satellites. At 
event M the receivers coincide in space-time but have different velocities. 

see Table I. The metric in Eq. Q is the solution of the linearized Einstein field equations 01- ^ n Eq. CJ, the 
coordinate time x° has no simple relation to the time kept by ideal clocks on the surface of the Earth. 

The coordinate time can be given a simple interpretation by transforming the metric in Eq. Q to rotating Earth- 
centered Earth- fixed (ECEF) coordinates y l , using the transformation 
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Note that the coordinate time in the rotating frame, y°, is equal to the coordinate time in the ECf frame x°. In these 
ECEF rotating coordinates, the metric is given by 
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From Eq. {3}, we see that there exist geopotential surfaces, (^(y 1 , y 2 , y 3 ) = c, where c is a constant, 



(4) 



t(y\y 2 ,y 3 ) = v-^ 2 ((y 1 ) 2 + ( y 2 ) 2 ) 



(5) 



where stationary clocks in the ECEF frame (that satisfy dy a = 0) have the same rate of proper time dr = ds/c with 
respect to coordinate time y° [l9j. In other words, ideal clocks located at the same value of geopotential 4>, have the 
same rate with respect to coordinate time y°. We have neglected the small cross-term 2VQ 2 R 2 / c 4 ~ 10 -21 . 
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Using the observation that clocks at a constant value of geopotential run at the same rate, it is advantageous to 
define 191 the new coordinate time t 



ct = X »=y»=(l + ^)f=ll + ^)x» (6) 



where (f> is the value of the geopotential <fi on the Earth's equator, at 6 — ir/2 and r = R: 

GM , 1 , 1 



For the values of the constants in Table I, the dimensionless magnitude of this term is (j> /c 2 — —6.96928 x 10~ 10 . 
Using the transformation in Eq. ©, the metric in Eq. Q becomes 



ds 2 = [1 + Jr - </>„)] (dy°) 2 + (l - ^) 2^dy 2 - y 2 dy')dy° 
+ (1 - ^) [(dy 1 ) 2 + (dy 2 ) 2 + (dy 3 ) 2 ] 



(8) 



Eq. (JHJ gives the space-time metric in ECEF rotating coordinates y l . Note that an ideal clock that is stationary in 
ECEF coordinates (with dy a = 0), has proper time 



dr = ds/c = — 
c 



1 + - 



1/2 



dy° 



(9) 



When this clock is located on the geoid, then 



0, and dr = dy° /c, so this ideal clock keeps coordinate time, 



x = y . Hence a good hardware clock that is on the geoid, and stationary with respect to the rotating Earth, can 
be used as a reference clock to keep coordinate time. Note that by Eq. © the coordinate time in rotating ECEF 
coordinates is the same as coordinate time in ECI coordinates, therefore, the same clock keeps coordinate time in the 
ECI frame, x°, and coordinate time in the ECEF frame, y°. 

Using the coordinate time transformation in Eq. JBJ, the metric in Eq. JQ) becomes 



ds = gijdx 1 dx 3 = — 



(dx ) 2 + (l- [(dx 1 ) 2 + (dx 2 ) 2 + (dx 3 ) 2 ] (10) 



Equation IjlUI) gives the metric in ECI coordinates. The coordinate time that enters into the metric, x°, is the time 
kept by ideal clocks on the geoid. This result was the goal of the time transformation given in Eq. 10. Note however, 
that in the ECI frame metric in Eq. (| 1 0|> . the proper time interval ds on a stationary clock in ECI coordinates (with 
dx a =0), is not equal to coordinate time interval dx° because in general V ^ 4> . The ECI coordinate metric, given 
in Eq. Q1U[I. is useful for computing the proper time dr = ds/c elapsed on-board a satellite clock, in terms of elapsed 
coordinate time. 



IV. RELATIVISTIC EFFECTS IN GPS 

A. Satellite Oscillator Frequency Offset 

The clocks in GPS satellites are at a higher gravitational potential than the clocks on Earth. As observed on the 
Earth, this difference in the gravitational potential causes the oscillators of GPS atomic clocks to appear to run fast, 
by fractional frequency 

Aa; = ^-^ 

LO C Z 

where the approximate gravitational potential (p = —GM/r, and r is the distance from the center of the Earth. Here, 
4>i and (f>2 are the potentials at the satellite and on the earth surface, respectively. For GPS, Aaj/ui ~ 5.28 x 10~ 10 . 
This effect causes the satellites clocks to run fast by 45 /is per day. This is often called a gravitational red shift, but 
actually, it is a blue shift (toward higher frequencies). 
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In addition to the gravitational frequency shift, the GPS satellites are moving. Consequently, as observed in the 
ECI frame, satellite oscillators exhibit time dilation due to their velocity v/c w 8.37 km/s. The time dilation effect in 
special relativity is given by 



A* = j At' 



(12) 



where At is the time interval in the ECI frame and At' is the proper time of a clock moving at speed v in the ECI 
frame. For GPS satellites, 7 » 1 — v 2 /(2c 2 ) w 8.33 x 10 _n . With respect to coordinate time in the ECI frame, the 
time dilation effect makes the satellite clocks appear to run slow by approximately 37 /us per day. 

The typical GPS atomic clock stability is 1 part in 10 -13 , so the effect of time dilation and gravitational red shift 
are each about 10 3 times larger, and therefore, both effects must must be taken into account. The net effect of time 
dilation and gravitational red shift is that the atomic clocks would run fast by 38 (is per day ( = 45 /lis - 7 £ts). This 
is a huge effect, which can be measured by the fact that 38 fis corresponds to a range error of 38,000 feet per day! 

The actual value of the combined effect of the gravitational potential and time dilation on the frequency of the 
satellite oscillator is computed using the metric in Eq. (|10fl . During a coordinate time dx°, the satellite moves a 
spatial distance dx a , a = 1,2,3. The proper time elapsed on the satellite clock, dr s — ds/c, is related to elapsed 
coordinate time, dx , by 



ds = c dr, = 



1 



bo) - (1 - ^f )vt 



1/2 



dx° 
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where V s is the Earth's gravitational potential, given in Eq. evaluated at the position of the satellite and v 2 = 
^ir^ftr is the square of the satellite velocity divided by c 2 . In Eq. Ijl3|l . we take the Earth's quadrupole potential 
to be zero, J2 = 0, which allows a circular orbit for the satellite and makes dr s /dx° independent of the polar angle 
8 of the satellite. Similarly, we approximate the satellite velocity v s by taking J2 = and assume a circular (zero 



eccentricity) orbit so that v 2 = GM/ac 2 
have 



where 



Expanding Eq. (|13|l to first order in small quantities, V s /c 2 and <fi /c 2 , we 
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(15) 



is the value of 4> when J2 = 0. Equation (|14fl gives an approximate expression for the rate of proper time on a 
GPS satellite in circular orbit, with respect to coordinate time x° in the metric in Eq. (|10fl . We have dropped terms 
V s v 2 /c 2 « 10~ 20 in Eq. JH|), since v 2 « 10~ 10 and V s /c 2 « 10~ 10 . 

The ratio of the frequencies of oscillators is inversely proportional to the clock rates, w COO rd/w s = dT s /dx°, where 
u cc ,ord is the number of cycles elapsed as counted per unit of coordinate time dx°, and u> s is the number of cycles 
elapsed as counted per unit of proper time dr s on-board the satellite. Therefore, the oscillators of clocks on-board 
the GPS satellites in circular orbit (taking J 2 = 0) have a frequency shift 0] 
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Using the value of JGM-2 constants |2C 
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If the clocks on-board GPS satellites were allowed to run freely, then Eq. I|14|) shows that during one day of elapsed 
coordinate time Ax°, a satellite clock would gain c(dT s /dx° — 1) x 24 Hours. Nominally, the GPS system is designed 
to transmit the digital P-code at 10.23 MHz. However, Eq. (|16[) shows that if the oscillator in the satellite were set to 
u! s /2ir — 10.23 MHz, then this code would appear to have a higher frequency, Lo COO rd = (1 + 5)oj s , as measured with 
respect to a clock that keeps coordinate time in the ECI frame. If a clock is stationary in the ECEF frame, and is 
located on the Earth's geoid, it keeps coordinate time x , and therefore the satellite clocks would appear to run fast 
to a GPS user on the geoid. Consequently, in the GPS, the technical specifications for the satellite clocks include a 
frequency "factory offset" that is applied prior to launch of the satellites. The actual (angular) frequency of P-code 
that is broadcast by the satellite clock is |2j| 



lo s = 27r(l - 5 ) fa 



(18) 



TABLE I: The numerical constants and values used in the calculations. Values taken from Ref. and |2(J. 



Symbol 


Definition 


Value 


Units 


GM 

R 

h 
a 

v s /c= y/GM/a/c 
GM/Rc 2 


Gravitational constant time Earth's Mass 

vflpinim m^ppd nf lip"Vit 

angular velocity of Earth rotation 
mean Earth radius at Equator 
Earth's quadrupole moment 
GPS orbit semimajor axis 
GSP satellite velocity / c 
gravitational potential scale 


3. 986004415 xlO 14 (JGM-2) 

2 QQ7Q2458x10 8 (exact definition 1 

7.2921151467xl0 -5 (WGS-84) 

6.3781363xl0 6 (JGM-2) 

1.0826269xl0 -3 (JGM-2) 

26561.75xl0 3 Ref. [19] 

1.29217xl0 -5 

6.95348xl0~ 10 


m7s J 
m/ s 

radian/s 

m 

1 

m 

1 

1 



where 6o — +4.46 x 10~ 10 and / = 10.23 MHz are fixed GPS constants. With this frequency correction applied, the 
satellite clocks approximately keep coordinate time x° in the ECI frame. 

However, satellite orbits cannot be made perfectly circular, so the GPS clocks are in a slightly eccentric orbit. This 
eccentricity of the orbit contributes an additional well-known effect on the time of the satellite clocks: the satellite 
clocks periodically speed up and slow down with respect to coordinate time x°. This effect is often called the "e Sin 
E effect" and depends on theposition of the satellite in its orbit. This clock correction is made in the GPS receiver 
in software, and is given by |22j 

2 

M r = —VGMae sinE (19) 

c 2 

where a is the semimajor axis of the satellite, e is the orbital eccentricity, and E is the eccentric anomaly along the 
orbit. For a typical upper limit of e = 0.01, At r w 23ns. 



B. Observed Frequency Shift 

Equation i|18|) gives the frequency offset that is applied to a GPS satellite clock/oscillator so that it approximately 
keeps coordinate time x°, with respect to a reference oscillator that is stationary in the ECEF frame and located on 
the geoid. Most users of GPS do not satisfy these two criteria. For example, a user that is stationary in the ECEF 
frame on the geoid is moving in the ECI frame, and a user in an aircraft is above the geoid. Therefore, at any given 
time, a user typically sees a frequency shift that is different from the "factory offset", given in Eq. (|18fl . The actual 
measured frequency shift of the satellite signal depends on user and satellite positions in the gravitational field of 
the Earth, and also on user and satellite velocities (not just user-satellite velocity differences). The reason that the 
frequency shift depends separately on user and satellite velocities is, of course, that space-time is not homogeneous 
(space-time is not flat) because of the Earth's gravitational field. 

Consider a satellite moving at velocity v s at gravitational potential V s . At event S, the satellite transmits an 
electromagnetic signal of proper frequency lo s , as measured with respect to a calibrated oscillator on-board the 
satellite. An observer at event O is at gravitational potential V and has a velocity v D . The observer measures a 
signal having frequency lu 0i which is different from lu s because he is in motion and at a different gravitational potential 
than the satellite. The emission event S and observer reception event O are connected by a null geodesic. The quantity 
lu /uj s — 1 is a 2-point scalar because it depends on two space-time points, S and O. A 2-point scalar transforms 
as a true scalar under separate transformations of coordinates at point S and at point O |23| . Using the metric in 
Eq. ljTU|) . a detailed calculation gives (see Appendix B) 

— - 1 = (v. - v ) • n [1 + ^ + \{v 2 - v 2 s ) + (v s • n) + (v s ■ fi) 2 - ttl - ft] 

+ ^ + Wo - V D + 7^7 [Mr ,r s ,v ,v s ) + ft(r s ,r ,v S) v )] (20) 

where r Q and r s are positions, and v c and v s are velocities (divided by c), of observer and satellite (source), respectively. 
In Eq. (|20|) . r = |r — r s |, v\ = v Q • v G , v 2 s = v s • v s , and V s and V D are the Earth's gravitational potential given in 
Eq. (J2J, evaluated at emission event 5* and reception event O, respectively. The unit vector fi connects the emission 
and reception events, S and O: 
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The dimensionless quantities a\ = GM/c 2 R and j3\ = Q 2 R 2 /2c 2 . The last term in the right hand side of Eq. lf2*U|) 
is given by the sum of two terms, which are related by interchanges of subscripts "s" and "o" . The scalar function 
h(r Q , r s , v Q , v s ) is given by 



h(r, 



1 



O ) 1 S •) v O ; 



1 



r s -v D (r„ • r s )(r c • v D ) 



(22) 



For a GPS user on the surface of the Earth, (V s — V )/c 2 w 5.3 x 10 10 , and the GPS satellite velocity (divided by c) 
is v s w 1.3 x 10 -5 . We also have that a x m 6.9 x KT 10 and ft » 1.2 x 10~ 12 . Therefore, in the derivation of Eq. QJJ, 
we have taken v s = 0(1), V s /c 2 = 0(2), V a /c 2 = 0(2), a x = 0(2), and ft = 0(2), where O(l) - 10" 5 , and I have 
dropped terms 0(4) — 1CT 20 . 

The measured frequency is due to three types of terms. First, there is a special relativistic Doppler effect that 
depends on the relative velocity of satellite and observer, contained in the terms proportional to (v s — v c ) • n, and 
the term \{v 2 — v 2 ), which comes from expanding the special relativistic 7-factors for satellite and observer. Next, 
there is a frequency shift due to the difference in gravitational potential of the observer and satellite, which is given 
by the stand alone term (V s — V )/c 2 . Finally, there are cross terms that depend on products of satellite and observer 
velocities and the Earth's mass M. 

As described in Eq. (JTSJl, the GPS satellites have a built-in "factory offset", <5 = 4.46 x 1CT 10 . From Eq. JSUJl, 
we see that the observed frequency shift of the satellite signal due to gravitational potential differences, (V s — V )/c 2 
is on the order of 10~ 10 and can vary depending on altitude of the observer. In addition to this frequency shift, 
there is a Doppler frequency shift that is much larger. For example, for a jet aircraft travelling along the equator 
at approximately 1000 km/hour with respect to the Earth's surface, the velocity (fraction of c) with respect to the 
ECI coordinates can be on the order v ~ 2 x 10~ 6 , which is a factor of 5 x 10 3 larger than the built-in "factory 
offset". Similarly, for an observer on-board a low-Earth orbit satellite at altitude 1000 km, whose orbit is in the 
plane of the equator, we have |v s — v G | ~ 10 -5 . This frequency shift is on the order of 10 4 times larger than the 
built-in "factory offset" . Of course the factory offset is applied to cancel out a secular effect-an effect that leads to a 
constantly increasing discrepancy in time between satellite and coordinate time clocks, while the Doppler (motional) 
effect has a more complicated time dependence. The point here is that the frequency of GPS signals, as seen by 
an observer, has huge frequency shifts due to observer motion. Below, we describe how these frequency shifts are 
essentially removed-so that the actual measurements made by a GPS receiver are independent of the velocity of the 
observer (GPS receiver). 



V. GPS SIGNALS: THE SPACE-TIME GRID 

The GPS satellites broadcast electromagnetic signals that set up a geometric space-time grid [24] . Users of the 
GPS that receive four satellite signals can uniquely identify their position in the space-time grid. This grid is created 
by discontinuities in the amplitude of the broadcast electromagnetic field. 

All the GPS satellites broadcast on two carrier frequencies in the L-band centered about: L\ « 1575.42 MHz 
and L2 « 1227.6 MHz. The GPS satellites have helical antennas that are constantly pointed toward the Earth 
center. Near the antenna axis (bore-site) the broadcast electromagnetic radiation field is right circularly polarized 
|15|. Superimposed on each satellite carrier signal is a unique code, or sequence of digital bits that identifies the 
satellite that is broadcasting. At a given time and spatial position in the comoving frame of the satellite, one vector 
component of the electric field can be written approximately as 

E = D(T*)s(T*)cos(quj s T* + <j>) (23) 

where lu s is the P-code transmission frequency given in Eq. I|18f) . t* is the hardware time kept on-board the satellite by 
its local clock, and cj) is a phase associated with phase noise in the signal. The integer q takes values q = 154 or q = 120, 
to produce the broadcast signal that is transmitted on the two carrier frequencies, L\ — 154 x oj s /2it ss 1575.42 MHz 
or L2 — 120 x oj s /2it « 1227.6 MHz. See Fig. [21 The L\ and L2 carrier frequencies are integer multiples of the code 
bit transmission rate, uj s /2tt. We distinguish between proper time r kept by an ideal clock, and hardware time, r*, 
which is kept by a clock that is a real physical device. The function s(t) takes the discrete values —1 and +1, and 
represents the digital P-code values and 1 in the pseudorandom code, which is unique to each satellite |2l|. There 
are exactly TV (as 10 12 ) values in the digital P-code sequence. The code sequence starts at midnight on Sunday and 
has a period of exactly one week: s(t + T) = s(t), where T is one week. Therefore, we can define a discrete phase $„ 
for the periodic function s(t) by 



*<») = 2 J^± 
N 



(24) 
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FIG. 2: Each component of the electric field, E(t*), broadcast by a GPS satellite has possible discontinuities in its amplitude 
at the times, t*T2 — Tjf where the pseudorandom P-code function s(r*) changes sign. These discontinuities propagate on 3-d 
hypersurfaces in space-time and define the space-time grid. 



where n = 1, 2, 3, N, is an integer that sequentially labels the bits in the code s(t). 

In Eq. 123|) . the function D(t*) is a digital navigation message that is broadcast at approximately 50 Hz. The 
message provides the satellite ephemeris and satellite clock corrections in the form of two coefficients, A and B. The 
clock corrections essentially provide the conversion from hardware time r* to coordinate time in the ECI frame, a; , in 
the form x° = A + Bt* . The functions s(r*) and D(t*) are timed so that changes between and 1 in s(r*) occur at 
the same point in time as in D(r*), i.e., the bit transitions (edges) in D(t*) align with those in s(r*). The function 
D(t*) provides the broadcast ephemeris of each satellite. 

In the satellite comoving frame, the pseudorandom code, s(r*), is broadcast at the factory adjusted angular fre- 
quency uj s given in Eq. I|18[i . as timed by the atomic clock on-board each satellite. The frequency of the oscillator 
(clock) on-board the satellite has been lowered from fo MHz to uj s /2it (see Eq. $181 ) so that, as seen from the 
ECI frame of reference, the discontinuities in the amplitude of the electromagnetic field occur at time intervals of 
approximately l//o with respect to coordinate time in the metric given in Eq. (|10f> 2^]. Since the code s(r*) is timed 
according to the satellite clock, the discontinuities in the code can occur at hardware times T*,n= 1, 2, 3, N, where 
n sequentially labels the (potential) code edge discontinuities since midnight Sunday, and iV is the number of code 
bits in the s(r*) code. All satellites have the same sequence of hardware times r* for possible code-edge-emission 
events, given by 

$(«) 2tt 

r: = — = -(n-l) (25) 
TV ^ 

where u> s is a constant frequency given in Eq. (|18fl and n — 1, 2, 3, • • • N is the bit number in the P-code sequence 
s(t*). In the ECI frame of reference, the coordinates of the world line of the antenna phase center for satellite s, 
xI(t*), can be parametrized by the satellite clock's hardware time, r*. The world line of a satellite, is approximately 
a geodesic, however, forces on a satellite, such as solar pressure, and antenna phase center position differing from the 
satellite center of mass, lead to an effective world line that differs from a geodesic. The world line of the antenna 
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FIG. 3: Discontinuities in the amplitude of the electromagnetic field at times, t*,t%, ■ ■ ■ t^, propagate on 3-d hypersurfaces in 
space-time that define the space-time grid. 

phase center is precisely determined by ground tracking stations. Future points on the world line are computed and 
expressed in terms of classical satellite ephemeris parameters and the A and B clock corrections. This information is 
uploaded to the satellite's navigation message, which is transmitted to GPS users through the digital sequence D(t*). 

It is well known that discontinuities in the emitted electromagnetic field tensor define invariant 3-d hypersurfaces 
called characteristics [I3,|2{|[27|. In terms of the world function of the space-time, for each satellite, these hypersurfaces 
are given by j2^| 

0(4^),^) =0 , n = 1,2, 3, •••AT (26) 



where x\ (r* ) is the world line of the satellite s parametrized by satellite hardware time r* , and the 3-d hypersurfaces 
are defined by their coordinates , j = 0, 1,2,3, that satisfy Eq. i|26|) and are on the forward light cone, so that 
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x° > £° (t* ) . Each hypersurface can be uniquely specified by t* , the hardware time of satellite "s" and the bit number 
n associated with the (discontinuity) new bit in the code s(t*). For flat space-time with a Minkowski metric, Eq. i|2(j|) 
reduces to 

i % (^-<«))(^~xi«)) = (27) 

where rjij is the Minkowski metric with nonzero diagonal components (—1, +1, +1, +1), and there is no sum on n. 
A detailed calculation of the world function entering in Eq. (|26|1 for metric in Eq. i|10|) gives 



^,4) = --(Ax ) 2 



2GM fl 2 R 2 



c 2 R 



where in Eq. (|28|l we use the following definitions, 



2GM 
c 2 |Ar| 



A( ri) r 2 )j +-(Ar) 2 

, 4 = (xP,,r 2 ), Ax° = 



2GM 
c 2 |Ar| 



A(ri,r 2 ) 



(28) 



x*2 and |Ar| 



ri -r 2 |, 



A(ri,r 2 ) = log 



tan(^) 



and 9\ and # 2 are defined by 



COS* 



(r 2 - ri 



r 2 - ri 



a = 1,2 



(29) 



(30) 



In Eq. H28|l . we have taken the small parameter J 2 = 0. 

Each GPS satellite broadcasts a set of 3-d hypersurfaces that form part of the coordinate grid, given in Eq. J2HJl. 
There are approximately 24 satellites in the GPS constellation, and all the hypersurfaces from these satellites comprise 
the GPS space-time grid that is used to label events in this space-time. The hypersurfaces are spaced approximately 
1/io ~ 97.75 ns in coordinate time and 29.31 m in space, see metric in Ea. l|10J) . Since the coordinate hypersurfaces 
are so far apart in space and time, an event in this space-time can be specified more accurately by interpolating the 
times t* at which discontinuities in electromagnetic field are actually emitted. For each satellite, s, an interpolated 
hardware time, t* and corresponding interpolated phase of pseudorandom code, $, can be associated with a real 
number r\ 



$ 2tt . 
— = — W' 



N 



1) 



(31) 



where cu s is given by Eq. fl%|) . r* < t* < r* +1 and n < r) < n+ 1, see Eq. (|2*5|) . The real number 77 is the interpolated 
value of the code bit integer number n. Associated with the interpolated hardware time, r*, there is a continuous 
family of coordinate 3-d hypersurfaces analogous to those in Eq. (|2*B|l : 



fi(x;(r„*),x») = 



(32) 



The continuous parameter r* labels the 3-d hypersurface that has code phase $ and is defined by coordinates x° that 
satisfy Eq. (|3*2*|) . On each 3-d hypersurface, the phase $ has the value 



$ = -J-T* 

N * 



(33) 



where uj s is given by Eq. (|18fl . The hardware time r* of an emission event at the satellite is related to coordinate 
time t s = a:° /c in the ECI frame of reference by the satellite clock correction At* 



t s = r* + At,; 

In terms of coordinate time of emission, the phase broadcast by satellite s is then given by 

$,(t,r) = ^( a: 2(t,r)- C Ar*) 



(34) 



(35) 



where lu s At* /^represents a phase offset due to the fact that satellite clocks keep hardware time, which is an approx- 
imation of coordinate time t s = x° s /c. The broadcast phase function satisfies the eikonal eauation|2f]j| 



dx l dxi 







(36) 
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where g l i are the contravariant components of the metric given in Eq. HlOfl (no sum on s). The phase function $ s 
depends on the invariant world line of the satellite, x\(t*). From Eci. (|36[) it is clear that under Lorentz transformations, 
or generalized coordinate transformations, the phase $ s (i, r) transforms as a scalar. The form of Eq.JSEJ shows that 
the wave vector associated with the phase & s , Ki — d& s /dx l , is a null vector. The wave vector Ki can be related 
to the covariant derivative of the world function of space-time, fi(T, i?), between emission event T = (tri^r) and 
reception event R = (t, r). The direction of the wavevector Ki is the same as the direction of the covariant derivative 
of the world function, f2i H , where the derivative is taken with respect to R. This can be seen from the identity 
f2 = g l; > (R)fli R Hj R = 0, since the geodesic connecting points T = (ir, i"t) and R — (t, r) is null. 

A world line is an invariant geometric quantity |2^|. For the limiting case of a satellite in Minkowski space-time, the 
time component t s of the satellite world line in Ea. l|35[) defines a scalar field, t s — t s (t, r), that is given by the implicit 
equation 

t a (t,r) =t- -|r-r,(t,(t,r))| (37) 
c 

where r s (i s ) is the ephemeris of the satellite. It is easy to check that this particular form for t s (t,r) in the phase 
function $ s (t, r) in Ea. (|35|l satisfies the eikonal Eq. l|36[l . 

For the actual case where we take into account the Earth's gravitational field, the light cone equation is given by 
fi(T, R) = 0, where f2 is the world function given in Ea. f28p . Writing the world function in Ea. (|28|l in the form 

Q = -1(1 + a)(x° - x°f + \{l + /3)(r - r s ) 2 (38) 

where a and (3 are small quantities, the scalar phase field $ s (r,t) in Ea. (|85fl can be written in terms of the function 
x® (x° , r) that is implicitly given by 



x° s (x°, F )=x°-(l+~(f3- a)) |r - r s (x° s (x°, r))| 



(39) 



where we have kept only linear terms in a and (3. Note that a and (3 are two-point functions that depend on 
space-time points T = (x®,r s ) and R = (x°,r). 

An event in this space-time can be uniquely labelled by four real hardware times, (rj* , t% , t£ , t\ ), or alternatively 
by dimensionless real numbers, (771, 772, 774), where each r\ s gives the continuous code parameter at emission time 
for satellite s. This system of coordinates has been studied by Synge who called them optical coordinates [23| . 
More recently, these same coordinates have been called GPS coordinates and their theoretical properties of have been 
explored in some detail[H[l3. 

At any time, a GPS user has more than four satellites in view, so in the real implementation of the GPS system of 
coordinates, they are multiple valued. Since the satellite signals are line-of-sight, a GPS user sees satellites rise and set 
on the horizon, and a different set of four satellites defines the 3-d hypersurfaces. There are currently approximately 
24 satellites in the GPS constellation, and 37 code sequence possibilities |2lJ, so at most, the numbers s — 1, 2, 3, • • • 37. 

In our discussion, we have neglected atmospheric effects. In practice, for users of GPS in the Earth's atmosphere, 
there exist significant propagation delays as well as frequency dispersive effects. The Earth's troposphere (atmosphere 
from the ground to approximately 10 km) causes the same time delay for both L\ and Li frequencies. On the other 
hand, the Earth's ionosphere (60 km to 700 km altitude) is dispersive at these frequencies, due to the presence of free 
electrons, and hence causes two different delays for the L\ and Li frequencies. A first principles treatment is possible, 
following Synge |23l |. However, in practice, simpler methods are used for correcting for both of these effects in GPS 
receivers |2l|. Atmospheric effects will not be considered further here. However, we note that signal propagation 
through the (frequency dispersive) atmosphere leads to a removal of the discontinuity in the electromagnetic tensor 
Fij. Hence, the dispersive effects of the atmosphere limit the accuracy of GPS for a ground user because the code 
edge is not sharp after signal propagation. The magnitude of this effect has not been investigated. Also, the finite 
size of the transmitting antenna also leads to a code edge that is not sharp [15| . 

VI. PSEUDORANGE MEASUREMENT 

A GPS receiver makes a special type of measurement, called a pseudorange measurement, in which the Doppler 
and gravitational frequency shifts are essentially removed (up to an additive constant) from the measurement. A GPS 
receiver internally replicates (generates) the satellite P-code s(t) for each satellite. The receiver determines which 
satellite signal it is receiving by matching the internally replicated P-code with the code received from the satellite. 
In what follows, we describe the pseudorange measurement made by a GPS receiver by using a mechanical analogue. 
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Our mechanical analogue of a GPS receiver has two code wheels, numbered 1 and 2, see Fig. 0] During operation of 
the receiver, both code wheels advance at the same angular velocity, Q R . The receiver will align the code incoming 
from the satellite with the code on wheel 2. The code arriving is like a bicycle chain that is to be placed onto the 
sprocket (code wheel 2). In order to align the codes, the receiver must match the rate (frequency) at which bits are 
arriving as well as match the duration of each bit. In general, the receiver has a non-zero velocity in the ECI frame 
and it is at a different gravitational potential than the transmitting satellite. At any instant in the receiver's comoving 
frame, the frequency of code bits received from the satellite is f R = lu r /2tt, where f R is the number of received bits 
per second of the code s(r*). To align the codes, this frequency lu r must be made equal to uj , which is the frequency 
seen in the comoving frame of the receiver, given in Eq.J30|). 

The temporal duration of the bits on code wheel 2 and the temporal duration of incoming bits from the satellite are 
made equal by adjusting the radius r w of code wheel 2. Code wheel 1 is adjusted to the same radius as code wheel 
2. If the code wheel has a circumference equal to 2irr w , the length per bit along the the circumference is 2wr w /N, 
where there are N bits in the code. In order to match incoming bits, the code wheel must have a linear velocity 
2irr w f R /N = r w Vt R . The angular velocity of the code wheel is then fl R = uir/N. This angular velocity of the code 
wheel, and consequently the frequency at which the replicated code is generated, u) R , will be adjusted continually and 
depends on the motion of the receiver and satellite, so that the replicated code is constantly aligned with the code 
arriving from the satellite: 

w R = lu (40) 

Therefore, a GPS receiver searches in time (aligns code phase on wheel 2 with incoming code) and in frequency 
(adjustment of r w to make temporal length of bits on wheel 2 equal to satellite bits). This alignment of the code is 
commonly referred to as "code tracking" or "code lock" . 

At any instant in time, code wheel 1 has a current value of phase, <& R , which is indicated by the pointer, see Fig. 
4. The value of the replicated phase <& R on code wheel 1 is a representation of the receiver's current hardware time, 
t r . The phase and receiver clock hardware time t r are related by 

*R = JfT* R (41) 

where u s is a system constant given by Eq. l|18|) . The system constant uj s is a conversion factor from hardware time, 
t r , to phase A key point is that during "code tracking" the phase $ R on code wheel 2 is adjusted to keep 
pace with the incoming code from the satellite, but this is done by advancing the time t Ri and not by changing the 
frequency u s , which is a fixed constant in the system. 

The replicated phase in the receiver can be related to the receiver's world line by writing the hardware time t r in 
terms of the coordinate time along the receiver's world line in the ECI frame, x R , 

4 = c(t r + At* r ) (42) 

where At r is the time correction to the receiver's hardware time, commonly called a clock bias. The receiver time, 
t r = T)T/N, is periodic, t r + T = t r , where T = 1 week (the period of the P-code s(t r )), and rj is a real number 
< f] < N, which corresponds to the interpolated value of the bit number n in the P-code s{t r ), see FigH This 
periodicity in t r does not introduce any ambiguity because it is assumed that a GPS user is located within one light 
week distance from the Earth. Using Ea. H42f> the receiver phase, can be related to the geometry by 

® R = ( X ° R ~ CAT ^ (43) 

On code wheel 2, the receiver has the replicated code (for a given satellite). This replicated code is rotated back 
and forth on code wheel 2 until it is aligned with the code phase & s (t R ,r R ) arriving from the satellite, where r R is 
the position of the receiver. As described above, at all times the phase on code wheel 2 is aligned with the incoming 
code phase, in both time (angle angle of the wheel) and frequency or bit duration (radius r w is adjusted). 

The pseudorange measurement made by a GPS receiver at space-time point (t R ,r R ) is the difference of the phases 
on code wheel 1 and 2 [3(j 

A$(t R ,r R ) = ^R-^s(tR,r R ) (44) 

where $> R is the curent value of the replicated code phase in the receiver (on code wheel 1) and $ s (t R , r R ) is the value 
of the broadcast phase Q(t,r) for satellite s (given in Ea. H35fl ) evaluated at the receiver at space-time measurement 
point (t R , r R ). The phase difference, A<fr(t Rl r R ), is a bonafide measurement because it is a comparison (a difference) 
between the value of the scalar field at the receiver, & s (t R ,r R ), and the replicated code phase, The measured 
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3>R 

FIG. 4: A mechanical analogue of a GPS receiver tracks the code broadcast from a GPS satellite. The receiver has two code 
wheels. The phase angle of code wheel 1 represents the current time in the receiver. The phase angle of code wheel 2 is kept 
aligned with the in-coming code sequence from the satellite. 
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difference, A&(t R ,r R ), is commonly expressed in units of length and is called a pseudorange measurement from 
receiver to satellite s: 

N 

p s = c—[$ R -$ s (t R ,r R )} (45) 

where p s is the measured pseudorange in units of length. 

The measured pseudorange, p s , is a scalar under generalized coordinate transformations (e.g., under Lorentz trans- 
formations) because the phase function & s (t, r) is a scalar field and $ R is an invariant that depends on the world line 
of the receiver. The pseudorange depends on both the world line of the satellite and the world line of the receiver, 
which are invariant quantities. The world line of satellite s, xI(t*) , enters into the definition of the scalar field & s (t, r), 
see Ea. (1351) . The world line of the receiver enters into the value of the phase of the replicated code, & R , see Ea. (l43|) . 
The world line of the receiver also enters into the evaluation of the scalar field $ s (i, r) at the position of the receiver at 
event R = {t R , v R ). Therefore, the pseudorange measurement, p s , is a two-point scalar that depends on the emission 
event T and the reception event R. The pseudorange is a particular case of quantities known as two-point tensors, 
which depend on two space-time points, and have tensorial transformation properties with respect to generalized 
coordinate transformations (that can differ) at each point |23j. Therefore, the pseudorange p s is a two-point scalar 
under separate generalized coordinate transformations at T = (ti^rr) and at R = (t R ,r R ) and can be labelled as 
p(tr, tt, t R , t_r);23]. The world function of space-time, fl(x\, x? 2 ) given in Eq.(|2EJ), is a well-known two-point scalar, 
and not surprisingly, the pseudorange and the world function have the same transformation properties. 

Two point-scalar fields enter into measurement situations whenever there is a field generated at space-time point 
T, such as an electromagnetic field, and a measurement or projection is done at the measurement event at space-time 
point R. In this sense, there is some similarity between a relativistic treatment of measurement processes and a 
quantum mechanical treatment, which is mentioned in the introduction. In the case of relativity, a measurement of 
a field is a projection of the field onto an observer's tetrad l2,Sl 1 3 1 1 - In quantum mechanics, traditionally projection 
operators are invoked in the definition of measurement [32 . l33l l34| . 



VII. TIME TRANSFER AND NAVIGATION IN SPACE-TIME 



The measured pseudorange, p s in Ea. H45|) . can be related to the geometric range between the emission event 
T = (tr, tt) and the reception event R = (t R , r R ) by using Ea.l(l3*]l for & R and evaluating the broadcast phase <E>(i, r) 
for satellite s (given in Eq.l |35(l 1 at the reception event R, leading to 



Ps = \t r -t t \ + cAt^ - cAt* r + A(r T , r B ) (46) 

where we have set Atj. = At*. In Ea. (|46[l . |i\r— tt\ is the geometric range between events T and R, At r is the clock 
correction to the receiver clock at event R, and cAr^ = x°(T)— cr*(T) is the satellite clock correction at event T. The 
satellite clock correction, Atj., is broadcast by the satellite in the navigation message contained in the function D(t*), 
see Eq. l|23ll . Conventionally, the pseudorange is modelled [9l ITol ITT| in flat space-time as p s = \r R — tt\ + cAt^~ cAt r . 
The term A(tt, tr) is a small correction due to the presence of the Earth's gravitational field that modifies space-time 
geometry near Earth. This correction depends on the mass of the Earth, M, the angular velocity of Earth rotation, 
O, and the Earth's equatorial radius, R, and is given by 

, 2GM ( .. , Irfl-r-rlX n 2 R 2 . 
A(r T ,r R ) = — — f A(r Tl r R ) - 1 H R 1 j - -^-|r fl -r T | (47) 

where the purely geometric function A(rx,r R ) is given in Ea. l29() . The quantity 2GM/c 2 w 0.89 cm is the gravita- 
tional radius of the Earth, and fl 2 R 2 \r R — yt\/c 2 w 4.8 x 10 -5 m, using |r^— tt\ ~ a — R, where a is the semimajor 
axis of the GPS satellites. The pseudorange in Eq, 146|l is based on the metric in the ECI frame of reference, which is 
given in Ea. (|10|) . 

If a user of the GPS knows their spatial coordinates in the ECI frame, r^, then a single pseudorange measurement 
p s to one satellite is sufficient to determine the users time t R from Eg. 146( 1 . In Eo. (|46|) yt = r s (tr), where r s (ty) is 
the broadcast satellite ephemeris evaluated at transmit time tr that satisfies 

Cl(tT,r s {t T ),t R ,r R ) =0 (48) 

where tt, t R , r R ) is the world function given in Eo. (|28|l . Therefore, Eos. (|46ll and l|48|) can be solved numerically 

for the user's reception event time t R , which is a coordinate time in the ECI metric in Ea. (|l(Jf) . 
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The more common case is that a user of the GPS wants to obtain time but he knows his spatial coordinates in the 
ECEF frame of reference, y = (y 1 ,y 2 ,y 3 )- The user's ECI frame coordinates, r^j = (x 1 , x 2 , x 3 ) , are given by 

r R = D(t T - 1 ) ■ y (49) 

where D(t — t$) is a time-dependent rotation matrix, which is equal to the unit matrix at epoch time t = to, see 
Eq.©. In this case, to obtain the user's time t R at the reception event, tr in Ea.l)4(i[l and l|48|l must be eliminated 
by use of Ea.l|4"9")l. Equations (^HJ), an d j£B) can be solved for tR. 

Navigation in space-time means the simultaneous determination of user position and time, i.e., determination of 
the space-time coordinates of the reception event R = (tn, r^). Navigation can be carried out by simultaneously 
measuring four pseudoranges, p s , to four different satellites, s = 1, 2, 3, 4. It is clear that four simultaneous equations 
of the form in Ea. H46|l can be solved for the four coordinates (tn, r^) that specify the user reception event in the ECI 
frame. Typically, a GPS user on Earth wants to know their coordinates in the ECEF frame of reference. The user 
must determine the coordinate time of emission for one satellite from Eq. (|48|l and use it in the transformation Eq. 1491) 
to determine user ECEF coordinates y = (y x ,y 2 ,y 3 ) |35j . 



VIII. RESOLUTION OF TWO RECEIVER EXPERIMENT 

If two receivers are at the same event M in space-time, and they are tracking the same four satellites, do the 
receiver's measure the same pseudorange? Also, do they obtain the correct space-time coordinates for the event Ml 
The resolution of the two receiver thought experiment, which was described in Section II, is now clear. The value 
of the phase received from each satellite, <& a (tR, tr), is the same for each receiver, independent of their velocities, 
because $ s (i, r) is a scalar field. However, the value of the hardware time for each receiver, and therefore the current 
replicated code phase, $r, is different for each receiver. This difference is due to the different clock bias Arjj of the 
two receivers. We can then say that each measured pseudorange is an invariant under coordinate transformations, 
however, the pseudoranges (to the same satellite) are different for each receiver because they depend on receiver world 
lines. For navigation purposes, even though the set of pseudoranges are different for each of the two receivers, the 
clock bias for each receiver is determined in the course of solving the four Eq. (|46|l . Therefore, while each receiver 
obtains a different set of four invariant pseudoranges to the same four satellites, each receiver obtains the correct 
space-time coordinates of the reception event M — (tR, Tr). 

IX. SUMMARY 

We used a two-receiver thought experiment in Section II to motivate the need to understand the transformation 
properties of the measured pseudorange, the quantity that is measured in the GPS. 

Starting from the weak field approximation for the metric of space-time in the vicinity of the Earth, we used the 
standard transformation to the rotating frame, and defined a new coordinate time so that coordinate time on the 
geoid surface is equal to proper time. Then we transform from the ECEF frame to the ECI frame, to obtain the metric 
that describes the ECI frame of reference near the Earth, given in Eq. IjlOl) . In this metric, due to a combination of 
gravitational and time dilation effects, the GPS satellite clocks appear to run fast with respect to coordinate time. 
We described the "factory offset" that is routinely applied to slow down the satellite clocks to (approximately) keep 
coordinate time in the ECI frame. In Section IV-B, for the metric in the ECI frame given in Eq.^UJ), we derived 
the apparent frequency of the satellite signal seen by an observer at an arbitrary position moving with an arbitrary 
velocity, due to Doppler and gravitational potential differences, see Eq.(|2UJ). 

We described the nature of the digital signal (P-code) that is broadcast by GPS satellites, and the space-time grid 
created by GPS satellites due to the discontinuities in the broadcast electromagnetic field. We computed the world 
function of space-time, given in Ea. (|28() . for the ECI metric in Ea. H10|) . Using this world function, we defined a scalar 
phase field 3> s (t, r), given in Eq. i|35|) . that satisfies the eikonal equation. This phase field is seen in the same way by all 
observers, independent of their state of motion. Using this phase field, and a mechanical analoque for a GPS receiver, 
we defined the measured pseudorange in Eq. l|45|l . and found that under generalized coordinate transformations it 
transforms as a two-point scalar, just like the world function. Within the geometrical optics approximation, there are 
no velocity effects on pseudorange measurements. Finally, we related the measured pseudorange to the geometry of 
space-time in Eq. <|4()|) . We obtained a small correction, given by A(rr,rij) in Ea.l|47|). to the conventional model of 
pseudorange. This correction is due to the curvature of space-time in the vicinity of the Earth. 
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APPENDIX A: CONVENTIONS AND NOTATION 



We use the convention that Roman indices, such as found on space-time coordinates x l , take the values i = 0, 1, 2, 3, 
and Greek indices take values a — 1,2,3. Summation is implied over the range of an index when the same index 
appears in a lower and upper position. If x l and x l + dx l are two events along the world line of an ideal clock, then 
the proper time interval between these events is dr = ds/c, where ds is given in terms of the space-time metric as 
ds 2 = —gij dx l dx J . We choose to have the signature +2, so when g^ is diagonalized at any given space-time point, 
the elements can take the form of the Minkowski metric given by 7700 = — lj VaP = &a/3- 



APPENDIX B: DOPPLER EFFECT IN EARTH'S GRAVITATIONAL FIELD 



A satellite at spatial position r s and travelling at velocity v s emits a signal in its comoving frame with frequency 
uj s . An observer at spatial position r Q and travelling at velocity v Q measures the frequency of this signal to be ui in 
his comoving frame. These frequencies are related by Ea. (|20|l . This relation can be derived by evaluating 23] 



UJ S 



n ls v l ° + n lo v i ° 



(Bl) 



where ili s and Q io are the covariant derivatives of the world function with respect satellite and receiver coordinates, 
x\ and x J , respectively, and V ts and V l ° are the components of the 4- velocities of satellite and receiver at space-time 
points S = x{ — (tT,rr) and O = x a — (tR,rR), respectively. The form of Ea. l|Bl|l together with the transformation 
properties of the world function, shows that the quantity uj /lj s is a two-point scalar under generalized coordinate 
transformations. We have computed the world function in Eq. (jBll) . ft = Sl(x\, xl), and it is given in Ea. l|28|l for the 
metric given in Ea. i|l(J|) . To obtain an explicit expression for uj /lo s we use an alternative derivation |23J. Consider 
the satellite and observer world lines, x % s {s a ) and x (s ), which are parametrized by proper times s s /c and s /c, 
respectively, see FigGO The lines AB and CD are null geodesies that connect satellite and receiver world lines, at 
successive cycles of the emitted signal. 

The frequency ratio lo /lo s is related to the proper times 



ds s 
ds 



(B2) 



where ds s is the proper time between AC and ds is the proper time between BD. Using the metric in Eci. (|10fl . the 
ratio of frequencies can be written as 



dfs 
ds n 



gooQg) +g a p(S)vfv : 
5oo(0) + 5q/3(0)wq 



1/2 



(B3) 



where the emission and reception events are connected by a null geodesic, so that £1(5, O) = 0- I 11 Eq. i|B3l) . g a /3(S) 
is the metric on the satellite world line at event 5 1 , g a p{0) is the metric on the observer world line at event O, and 
i>™ and Vg are the velocity components of satellite and receiver, for a = 1,2,3. For example, the satellite velocity 
is vf — Aa;"/Ax^, where Ax® is the coordinate time between events AC on the satellite world line. In the limit 
Ax® — ► (and therefore Ax® — > 0), using the world function in Eq.J2HJ) to connect points AC and BD, a lengthy 
calculation of Ea. (|B3|) using the metric in Ea. ljlT)|) leads to Ea. lf2*U|l . The term Ax® /Ax® in Ea. (|B3|) leads to linear 
velocity terms in Ea. l(2T)|l . 
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